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Military  Construction,  Office,  Chief  of  Engineers,  U.  S.  Army,  under 
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Unit  005,  Study  of  Engineering  Classification  of  Cohesionless  Soils. 
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Mr.  C.  L.  McAnear,  Chief,  SMD. 
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I.  INTRODUCTION 


The  purpose  of  this  study  was  to  determine  fabric  characteristics  of 
undisturbed  sand  samples  from  Niigata,  Japan.  The  results  of  these  studies 
in  combination  with  results  of  completed  and  continuing  studies  on  the 
sampling  and  testing  of  sands  from  Niigata  (Silver  and  Ishihara,  1977) , as 
shown  in  Table  1,  have  as  their  objective  development  of  improved  under- 
standing of  the  liquefaction  behavior  of  sands  in-situ.  It  has  already  been 
established  that  liquefaction  properties  of  sand  samples  compacted  in  the 
laboratory  bear  some  relationship  to  method  of  sample  preparation  (Ladd, 

1974;  Mulilis  et  al,  1977)  and  to  the  sand  fabric  as  characterized  by  inter- 
particle contact  orientations  and  particle  long  axis  orientations  (Mitchell, 
1976).  Previous  work  by  Oda  (1972,  1976),  Mahmood  (1973)  and  Mitchell  (1976) 
has  shown  that  the  static  strength  and  compressibility  of  sands  and  silts 
also  depend  upon  soil  fabric.  It  is  therefore  suspected  that  in-situ 
liquefaction  behavior  of  sand  deposits  in  response  to  seismic  loading  may 
be  related  to  some  degree  with  the  sand  fabric. 

A fabric  analysis  was  performed  on  undisturbed  sand  samples  from  two 
sites  in  Niigata,  Japan,  as  shown  in  Fig.  1.  One  site,  the  "River  Site", 
experienced  significant  distress  from  liquefaction  during  the  1964  earthquake, 
while  the  other,  the  "Road  Site,"  showed  no  surface  evidence  of  liquefaction 
during  the  earthquake.  Site  profiles  are  presented  in  Fig.  2.  A complete 
description  of  the  site  profiles  and  topography  has  been  presented  by  Silver 
and  Ishihara  (1977) , who  state  that  it  is  generally  believed  that  liquefaction 
at  the  River  site  only  occurred  in  the  upper  9 meters  (30  ft  ) of  the  soil 
profile. 
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A large  number  of  samples  were  obtained  at  the  Niigata  site  using 
both  an  Osterberg  Sampler  and  a newly  developed  Japanese  Large  Diameter 
Sampler.  Only  Osterberg  samples  were  chosen  for  use  in  this  study  as  the 
Japanese  sampling  method  was  not  used  at  the  Road  site.  Depths  and 
identification  numbers  are  indicated  in  Fig.  2 for  the  eleven  samples  chosen. 

It  is  important  to  emphasize  at  the  outset  that  as  the  sand  at  the 
River  site  liquefied  during  the  1964  earthquake,  present  fabric  differences 
between  samples  from  the  Road  and  River  sites  cannot  be  related  to  lique- 
faction potential  prior  to  that  earthquake.  They  may,  however,  relate  to 
present  liquefaction  potential. 
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II.  SAMPLING  PROCEDURE 

Introduction 

[ 

An  important  aspect  of  this  overall  project  was  the  development  of 
new  soil  sampling,  sample  handling  and  sample  storage  techniques  in  an  effort 
to  minimize  sample  disturbance  between  the  field  and  the  laboratory.  A 
complete  description  of  these  procedures  is  provided  by  Silver  and  Ishihara 
(1977).  In  the  following  section,  details  of  sampling  technique  will  be 
discussed  only  as  they  pertain  to  the  subsequent  fabric  study. 

Sampler  Description  and  Operation 

The  samples  used  in  this  study  were  obtained  by  means  of  a U.S. 

Osterberg  Sampler,  a hydraulically  operated  piston  sampler  (Osterberg,  1952) . 

A schematic  illustration  showing  the  operation  of  the  Osterberg  Sampler  is 
presented  in  Fig.  3. 

Sampling  at  Niigata  was  performed  in  a mud  filled  borehole  100mm  (4  in.) 
in  diameter.  The  hole  was  advanced  with  a fish- tailed  bit,  which  was  modi- 
fied with  baffles  that  directed  the  drilling  fluid  away  from  the  bottom  of 
the  borehole  to  reduce  disturbance.  Fresh  drilling  mud  was  introduced  into 
the  borehole  as  the  sampler  was  withdrawn  to  prevent  caving.  When  the  sampl- 
ing tube  cleared  the  riser  casing,  a small  amount  of  soil  was  removed  from 
the  core  tube,  an  expandable  0-ring  packer  installed  in  the  tube  and  locked 
to  prevent  the  loss  of  sand.  A small  hole  was  then  drilled  below  the  piston 
to  release  the  vacuum  developed  by  the  sampler  and  to  make  it  easier  to 
remove  the  sample  tube.  After  removing  the  sample  tube,  the  dimensions  of 
the  sample  were  recorded. 

The  bottom  packer  was  equipped  with  small  vertical  drainage  holes. 

When  the  vacuum  under  the  piston  was  released,  water  from  the  soil  voids 
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could  flow  out  of  the  sample.  Samples  were  drained  for  24  hours  in  a 
vertical  position  to  prepare  them  for  subsequent  freezing. 

The  horizontal  orientation  of  the  samples  was  not  recorded  on  the 
sampling  tubes.  This  means  that  a truly  three  dimensional  comparison  of 
sample  fabrics  from  different  locations  cannot  be  made. 

Specimen  Handling 

After  24  hours  of  free  drainage,  the  samples  were  quickly  frozen  in 
liquid  nitrogen  in  the  field  and  stored  in  dry  ice  while  they  were  trans- 
ported from  Japan  to  the  laboratory  at  the  University  of  California.  In  the 
laboratory,  the  samples  were  stored  in  a commercial  freezer  at  -10°C. 

The  quality  of  the  frozen  specimens  was  checked  by  determining  the 
specimen  volume  before  and  after  freezing  (Silver  and  Ishihara,  1977) . No 
measurable  difference  was  noted,  although  an  average  one  percent  volume 
decrease  accompanied  subsequent  thawing.  The  fact  that  freezing  did  not 
change  the  dimensions  of  the  sample  led  Silver  and  Ishihara  (1977)  to 
conclude  that  the  drainage  procedure  effectively  cleared  the  soil  voids  of 
excess  water  and  only  provided  enough  water  at  the  intergranular  contacts 
to  provide  particle  binding  when  the  specimen  was  frozen. 

The  lack  of  measurable  sample  deformation  during  freezing  indicates 
that  in-situ  void  ratio  was  preserved,  but  since  sands  at  equivalent  void 
ratios  may  possess  radically  different  fabrics,  the  possibility  of  an 
induced  change  in  fabric  as  a result  of  freezing  does  exist.  To  investigate 
this  possibility  further,  a limited  study  was  made  by  Guzikowski*  using 
Monterey  No.  0 sand.  It  was  found  that  for  samples  prepared  by  moist  tamping 

*Guzikowski,  F.  J. , Unpublished  Research  Report,  Department  of  Civil 
Engineering,  University  of  California,  Berkeley,  1977. 
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freezing  had  no  statistically  significant  effect  on  the  particle  long  axis 
orientations  for  relative  densities  of  50  and  80  percent.  The  distribu- 
tions of  interparticle  contact  normal  orientations  were  not  statistically 
different  for  frozen  and  unfrozen  samples  at  a relative  density  of  80  per- 
cent. For  samples  at  a relative  density  of  50  percent,  however,  it  was 
found  that  on  the  order  of  15  percent  of  the  interparticle  contacts  were 
disturbed,  as  evidenced  by  changes  in  the  interparticle  contact  distribu- 
tion function. 

Soil  water  suction  - water  content  determinations  were  also  made  by 
Guzikowski  on  samples  that  had  been  previously  frozen  for  comparison  with 
results  for  unfrozen  samples.  From  these  data  estimates  of  pore  size 
distribution  could  be  made.  No  significant  difference  was  observed  for 
samples  compacted  to  80  percent  relative  density.  Small  differences, 
attributable  to  either  real  pore  size  differences  or  to  sample  seating  in 
the  test  cell,  were  observed  for  samples  at  50  percent  relative  density. 

It  seems,  therefore,  that  freezing  may  cause  fabric  changes  for 
samples  at  low  relative  densities,  although  preliminary  results  by  Walberg 
(1977)  indicate  that  freezing  may  not  influence  cyclic  triaxial  results. 

It  should  be  noted  further  that  as  the  sand  at  Niigata  was  sampled 
prior  to  freezing  rather  than  frozen  in-situ  prior  to  sampling  (coring  would 
be  required) , the  observed  fabric  may  differ  from  the  actual  fabric  in  the 
ground. 
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III.  SAMPLE  PREPARATION 


Introduction 

The  frozen  Osterberg  samples  as  received  at  the  laboratory  were 
2.8  inches  (71  ran  ) in  diameter  and  approximately  7 inches  (178  mm  ) in 
length.  The  top  portion  of  each  sample  was  used  for  the  preparation  of 
thin-sections,  as  well  as  sections  for  X radiographs,  while  the  bottom 
portion  was  used  in  measuring  the  relationship  between  water  content  and 
soil-water  suction.  Thin-section  and  X radiograph  section  preparation 
necessitated  resin  impregnation  of  dried  samples,  while  the  dimensions  of 
the  "Tempe-Cell"  soil  suction  device  dictated  that  the  samples  be  cut  and 
trimmed  prior  to  testing.  In  an  effort  to  preserve  the  fabric,  samples 
were  trimmed  while  frozen. 


Sample  Cutting 

To  obtain  portions  of  the  samples  for  both  resin  impregnation  and 
frozen  trimming,  the  san^sle  tubes  were  cut  with  a pipe  cutter  into  two 
sections  of  approximately  5 and  2 inches  (127  and  50  mm  ) in  length.  Once 
the  tubes  were  cut,  a fine  toothed  hacksaw  was  used  to  cut  the  sample  itself. 
The  entire  cutting  operation  could  be  performed  in  a few  minutes.  In  each 
case,  visual  inspection  indicated  that  no  significant  thawing  of  the  sample 
had  occurred  outside  the  immediate  cutting  zone. 


Sample  Thawing 

The  portion  of  the  sample  to  be  impregnated  was  fitted  with  a 
specially  fabricated  base  cap  equipped  with  a threaded  drainage  port  and 
wire  mesh  filter.  The  sample  tube  was  placed  vertically,  and  the  sample 
height  was  measured  at  different  points  to  the  nearest  .01  in.  (0.25  mm  ). 
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The  sample  was  allowed  to  thaw  and  drain  freely  overnight.  No  measurable 


change  in  sample  height  was  observed  as  a result  of  thawing.  It  was  noted, 
however,  that  on  the  order  of  0. 5-1.0  ml.  of  water  had  flowed  out  of  the 
drainage  port  during  thawing.  This  outflow  may  be  traced  to  two  possible 
sources:  1)  atmospheric  condensation  on  the  cold  sample  and  sample  tube, 

2)  additional  drainage  from  the  sample  which  could  be  attributable  to 
incomplete  drainage  before  freezing,  to  a change  in  the  pore  size  distribu- 
tion during  freezing,  handling  or  thawing,  and/or  to  a reduction  in  the 
water  adsorption  tendencies  of  the  sand  grains. 


Resin  Impregnation 

After  thawing,  the  moist  samples  were  dried  in  a forced  air  oven  at 
60°C.  Warm  air  was  drawn  slowly  through  the  sample  by  connecting  the  base 
cap  to  a vacuum  of  about  1 psi  and  allowing  warm  air  to  enter  through  the 
top.  The  samples  were  allowed  to  dry  in  this  manner  for  8 hours  or  more. 

The  sample  to  be  resin  impregnated  was  then  cooled,  fitted  with  a 
plexiglas  cap  and  subjected  to  a vacuum  of  5 psi.  Koppers  polyester  resin 
No.  1086-4,  mixed  with  24  percent  (by  volume)  styrene  (Koppers  No.  VM-901) 
to  reduce  viscosity  and  1 percent  (by  volume)  organic  peroxide  (Witco 
Chemical  HI-Point-180)  as  catalyst,  was  introduced  at  the  bottom  of  the 
sample.  After  saturation  the  sample  was  left  undisturbed  until  initial  set 
of  the  resin  had  occurred.  The  sample  was  then  placed  in  a forced  air 
oven  at  60°C  for  24  hours  to  permit  the  resin  to  cure. 

After  curing,  the  samples  were  removed  from  the  oven  and  allowed  to 
cool.  The  sample  tubes  were  cut  longitudinally  with  a bandsaw  and  the 
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Thin  Sections 


Three  thin  sections  were  cut  from  the  central  portions  of  the  impreg- 
nated samples  at  orthogonal  orientations  as  shown  in  Fig.  4.  By  determining 
grain  and  contact  orientations  in  each  of  these  thin  sections,  their  dis- 
tributions in  all  directions  within  the  sample  could  be  determined. 

Petrographic  thin  sections,  about  60  ym  thick,  20  mm  wide,  and  30  mm 
long,  were  prepared  using  standard  methods. 

X-Radiograph  Sections 

The  portions  of  the  resin  impregnated  samples  remaining  after  thin 
section  preparation,  were  used  for  X-radiograph  inspection.  The  largest 
remaining  pieces  of  each  sample  were  selected,  and  horizontal  as  well  as 
vertical  sections  about  0.2  inch  (5  mm)  thick  were  cut  by  means  of  a 
circular,  water  cooled,  rock  saw.  Sawing  was  completed  in  stages  to  pre- 
vent melting  of  the  resin.  As  horizontal  orientations  of  the  samples  were 
not  determined,  the  orientations  of  the  vertical  sections  were  selected  to 
provide  maximum  dimensions  for  each  section. 


Trimming  of  Frozen  Samples 

The  portion  of  the  sample  to  be  used  in  the  pore  size  distribution 
study  was  trimmed,  while  frozen,  to  the  required  dimensions.  The  sample 
tube  was  cut  lengthwise  with  a band  saw.  The  heat  generated  in  this  proce- 
dure produced  a thawed  zone  in  the  vicinity  of  the  cut.  The  sample  was 
returned  to  the  freezer  and  care  was  taken  to  note  the  extent  of  thawing  to 
avoid  the  subsequent  use  of  disturbed  portions  of  the  sample. 

The  frozen  sample  was  trimmed  to  the  dimensions  indicated  in  Fig.  5 
with  the  aid  of  teflon  coated  aluminum  templates.  A hacksaw  blade  was  used 
to  cut  the  sample  into  a disk  and,  with  the  templates  as  a guide,  metal 
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files  were  used  to  trim  the  sample  to  size.  The  trimming  procedure  was 
performed  in  short  stages  and  the  sample  was  returned  to  the  freezer  for 
at  least  one  hour  between  each  stage. 

After  trimming,  the  dimensions  and  weights  of  the  frozen  samples 
were  recorded. 
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IV.  SAND  PROPERTIES 

Introduction 

The  Niigata  sand  is  a fine  to  medium  sand  composed  primarily  of 
angular  to  subangular  quartz  particles.  The  Shinano  River,  Fig.  1,  has 
meandered  over  the  site,  with  alternate  cycles  of  scour  and  deposition  of 
silt  and  sand.  At  both  the  River  and  Road  sites,  the  sand  at  depths 
greater  than  4 m consists  of  fluvial  river  deposits.  The  sand  from  1.0  to 
4.0  m depth  at  the  River  site  and  1.5  to  4.0  m depth  at  the  Road  site 
was  dumped  through  water  and  not  subsequently  densified.  Uncompacted  fill 
was  placed  above  to  the  present  ground  surface. 

The  sample  properties  which  were  quantitatively  examined  in  this 
fabric  study  included  void  ratio,  grain  size  distribution  and  particle  axial 
ratio  distribution.  Portions  of  each  sample  which  were  disturbed  or  thawed 
during  cutting  and  trimming  were  retained  for  grain  size  and  axial  ratio 
measurements.  The  weights  and  dimensions  of  thawed,  dried  samples  were  used 
in  conjunction  with  specific  gravity  values  obtained  by  Silver  and  Ishihara 
(1977)  to  calculate  void  ratios. 

Grain  Size  Distributions 

Grain  size  distributions  were  obtained  by  dry  sieving  the  oven-dried 
portions  of  the  sample  retained  after  cutting  and  trimming.  Grain  size 
distributions  of  the  River  site  samples  are  presented  in  Fig.  6.  Grain  size 
distributions  of  the  Road  site  samples  appear  in  Fig.  7.  The  samples  from 
both  sites  are  seen  to  fall  into  the  medium  to  fine  sand  size  range. 

The  mean  particle  size,  D_n,  and  the  uniformity  coefficient,  D /D  , 
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are  plotted  for  each  site  as  functions  of  depth  in  Fig.  8.  Results  obtained 
by  Silver  and  Ishihara  (1977)  are  also  plotted  for  comparison. 
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The  mean  particle  size  of  the  River  site  samples  exhibits  a uniform 


increase  with  depth.  Values  of  Dc  obtained  from  Road  site  samples  are 


slightly  lower  and  more  variable  than  those  of  the  River  site  samples.  The 


mean  particle  sizes  measured  by  Silver  and  Ishihara  (1977)  are  seen  to 


agree  closely  with  the  values  obtained  in  the  present  study.  The  uniformity 


coefficients  of  the  Road  site  samples  are  slightly  lower  than  those  of  the 


particle  size  and  coefficient  of  uniformity  for  the  sand  at  the  Road  site 


would  suggest  a greater  susceptibility  to  liquefaction.  Even  for  a given 


input  ground  motion,  all  things  would  not  be  equal,  however,  as  the  water 


table  is  almost  2 m deeper  at  the  Road  site  than  at  the  River  site 


The  uniformity  coefficients  calculated  by  Silver  and  Ishihara  (1977) 


Fig.  8,  are  consistently  lower  than  the  values  measured  in  this  study.  It 


should  be  noted  that  the  grain  size  distributions  presented  by  them  were 


determined  using  Japanese  standards.  This  discrepancy  in  measured  uni 


Axial  ratio,  AR,  as  discussed  by  Oda  (1976)  is  an  index  of  particle 


shape  which  is  calculated  by  dividing  the  shortest  dimension  (L  ) of  a 


projected  particle  image  by  the  length  (L1)  of  a rectangle  tangent  to  the 


particle  image  outline,  as  illustrated  in  Fig.  9.  Thus,  axial  ratio  is 


seen  to  be  the  inverse  of  the  particle  length/width  ratio,  L/W,  employed  in 


earlier  studies  by  Mahmood  (1973)  and  Mitchell  et  al.  (1976) . Axial  ratio 


was  selected  as  a particle  shape  index  for  use  in  this  study,  as  the  value 


which  AR  can  assume  is  limited  to  the  range  0 < AR  <.  1.  When  the  axial  ratio 


approaches  unity,  the  projected  particle  image  is  close  to  circular.  The 


tic  of  particle  shape  with  similar  meaning  to  particle  sphericity  as 


discussed  by  Oda  (1976) 


metic  mean  of  a number  (approximately  200)  of  axial  ratio  observations.  An 


investigation  by  Oda  and  Koishikawa  (1977)  indicates  that  statistically 


which  exhibit  a value  of  AR  less  than  approximately  0.65 


than  0.65,  indicating  that  a statistically  significant  study  of  particle 


exhibit  a slightly  wider  range  of  particle  axial  ratio  values  and  possess  a 


higher  average  value  (.62  vs  .59)  than  the  Road  site  samples 


Sample  Void  Ratio 


Dry  unit  weight  values  were  determined  by  measuring  the  dimensions  of 


the  thawed  samples  which  were  prepared  for  resin  impregnation.  The  dried 


sample  was  weighed  while  still  in  the  sample  tube,  and  the  weight  of  the 


length)  were  subtracted.  The  accuracy  of  calculating  tube  weight  from 


measured  tube  length  was  determined  by  applying  the  procedure  to  a number  of 


scrap  sample  tube  lengths.  It  was  found  that  the  weight  of  a five  inch 


(127  mm)  length  of  tubing  could  be  predicted  to  within  .01  lb.  (4.54  gm) 


less  than  two  percent  in  the  computed  void  ratio 


Void  ratios  calculated  in  this  manner  are  presented  in  Fig.  12 


together  with  additional  void  ratios  measured  during  the  preparation  of 
samples  for  soil  suction  testing.  Results  obtained  by  Silver  and  Ishihara 
(1977)  are  also  shown.  The  void  ratios  of  samples  used  in  this  study  are 
seen  to  be  consistent  with  the  range  of  values  obtained  by  Silver  and 
Ishihara  from  Osterberg  and  Large  Diameter  samples  at  similar  depths.  The 
data  indicate  that  the  in-situ  void  ratios  for  the  River  and  Road  sites  are 
distributed  erratically  with  depth.  Void  ratios  at  the  Road  site  are  seen 
to  be  generally  higher  than  void  ratios  at  the  River  site. 

Insufficient  sample  quantities  were  available  to  determine  maximum 

and  minimum  void  ratios  according  to  ASTM  standard  D2049.  Limiting  void 

ratios  (e  , e . ) as  presented  by  Silver  and  Ishihara  (1977)  were  there- 
max  mm 

fore  used  to  calculate  the  relative  density  of  the  samples  studied.  The 
calculated  values  are  subject  to  a degree  of  uncertainty  since  the  Japanese 
procedure  used  to  measure  efflin  produces  samples  with  a higher  maximum  density 
than  is  attainable  using  the  more  usual  ASTM  procedure.  An  additional 
uncertainty  in  the  calculated  relative  density  values  arises  from  the  fact 
that  natural  sample  variation  may  cause  the  limiting  void  ratios  as  deter- 
mined by  Silver  and  Ishihara  to  differ  from  the  limiting  void  ratios 
appropriate  for  the  fabric  study  samples. 

Relative  densities  for  the  two  profiles  are  presented  as  functions  of 
depth  in  Fig.  13.  Results  obtained  by  Silver  and  Ishihara  (1977)  from 
Osterberg  and  Large  Diameter  samples  are  also  shown.  The  values  obtained 
from  the  fabric  study  samples  are  seen  to  be  consistent  with  the  range  of 
values  obtained  by  Silver  and  Ishihara.  The  data  indicate  that  the  in-situ 
relative  density  values  are  low  at  both  sites,  typically  ranging  from  20  to 
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60  percent.*  Relative  densities  of  River  site  samples  measured  in  this 
study  are  seen  to  be  slightly  larger  than  those  of  the  Road  site  samples. 
However,  the  wide  range  of  measured  by  Silver  and  Ishihara  indicates 
that  sample  density  is  too  variable  to  warrant  the  drawing  of  definite 
conclusions  from  the  values  obtained  from  the  samples  used  in  this  study. 


*Were  the  ASTM  standard  procedure  used  to  determine  the  maximum  density,  a 
larger  value  of  e^n  would  have  been  obtained,  resulting  in  somewhat 
greater  values  of  relative  density. 


Fabric  elements  studied  were  the  three-dimensional  distributions  of 
normals  to  interparticle  contact  planes,  the  orientations  of  apparent 


particle  long  axes  in  the  plane  of  the  thin  section,  and  the  distribution 
of  pore  sizes  as  inferred  from  the  measurement  of  soil  suction-water  content 
relationships.  Other  fabric  elements  such  as  void  variability  within 
samples,  particle  clustering  and  particle  size  segregation  were  not  examined 
quantitatively. 

Universal  Stage  Microscope 

A microscope  equipped  with  a universal  stage  was  used  to  determine 
particle  and  contact  orientations.  The  universal  stage,  shown  schematically 
in  Fig.  14,  allows  independent  rotation  of  the  thin  section  about  three 
orthogonal  reference  axes.  The  thin  section  is  coated  with  an  oil  possess- 
ing  the  appropriate  index  of  refraction,  and  covered  by  the  upper  glass  hemis- 
phere (H  in  Fig.  14)  set  in  the  metal  mount,  M.  The  thin  section  is  aligned 
and  moved  with  the  aid  of  the  sledge  R which  moves  horizontally  in  a slot 
cut  in  the  hemisphere  mount.  Rotation  can  be  made  about  the  NS,  E IK,  and 
vertical  axes,  and  the  amount  of  rotation  is  measured  on  the  appropriate 
vernier. 

The  use  of  the  universal  stage  microscope  is  described  in  detail  by 
Turner  and  Weiss  (1963) . 


Sampling  Methods 

Selection  of  particles  for  the  study  of  long  axis  orientations  and 
contact  plane  orientations  in  a given  thin  section  can  be  considered  as  a 
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two-dimensional  sampling  problem,  because  the  thickness  of  the  thin  section 
is  very  small  compared  to  its  length  and  width.  Several  methods  of  sampling 
two-dimensional  populations  are  available  (Kel lerhals  and  Bray,  1971;  Koch 
and  Link,  1970;  Miller  and  Kahn,  1962,  Pincus,  1953): 

1.  Sample  the  entire  population. 

This  is  the  only  method  that  can  describe  the  true  orientation  of  the 
entire  population.  This  method  is  impractical  for  the  typically  large 
populations  of  the  thin  sections  used  in  this  study. 

2.  Grid  sampling. 

A grid  is  established  over  the  sample  and  only  those  grains  or  con- 
tacts immediately  below  the  grid  points  are  counted.  A variation  of  this 
method  is  to  assign  a number  to  each  grid  block  and  to  count  the  entire 
populations  within  randomly  selected  grid  blocks. 

3.  Transect  sampling. 

Every  member  of  the  population  that  falls  under  a straight  line  or 
narrow  strip  across  the  sample  is  counted. 

4.  Pseudo-random  sampling. 

Arbitrarily  chosen  members  of  the  parent  population  are  measured. 

This  is  most  susceptible  to  sampling  bias. 

The  transect  sampling  method  was  used  in  this  investigation,  because 
(1)  the  population  was  too  large  to  sample  in  its  entirety,  (2)  no  grid 
attachment  was  available  for  the  microscope  and  (3)  it  was  necessary  to 
avoid  sampling  bias.  Traverses  about  eight  or  nine  grains  wide  and  about 
one  inch  long  were  chosen  parallel  to  the  long  axis  of  the  thin  sections. 
Every  grain  or  contact  within  these  zones  was  counted  until  a statistically 
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significant  number  of  observations  had  been  recorded  as  discussed  in  more 


Apparent  Longest  Axis  Orientations 


The  orientation  of  a particle  is  represented  by  the  inclination  of 


its  true  long  or  short  axis  with  respect  to  fixed  reference  axes.  In 


practice  the  true  orientation  of  small  particles  is  estimated  by  analyzing 


the  orientation  of  apparent  longest  axes  as  seen  in  thin  sections.  The 


three-dimensional  orientations  of  grains  can  estimated  from  the  orienta 


mined  in  two  vertical  and  one  horizontal  section  for  each  sample 


Orientations  were  determined  for  about  200  grains  in  each  thin  section 


to  the  vertical  direction.  For  the  horizontal  thin  sections 


with  0=0 


Statistical  analyses  of  the  long  axis  orientation  data  were  performed 


using  a vector  method  and  test  of  significance  presented  by  Curry  (1956) 


According  to  this  method  each  long  axis  observation  is  considered  to  be  a 


vector  of  unit  magnitude  in  the  measured  direction.  No  distinction  is  made 


A resultant  vector  which  is  calculated  from  these  unit  observation 


vectors  may  not  reflect  the  true  central  tendency  of  the  distribution.  If 


the  180°  range  lies  in  an  eastern  semi-circle,  the  distribution  has  no 
west  components  at  all.  North  components  would  tend  to  cancel  south 
components  so  the  resultant  vector  would  always  possess  a strong  easterly 
tendency  even  if  the  true  central  tendency  lay  close  to  the  north-south 
direction.  To  remedy  this  bias  the  angles  of  the  observation  vectors  are 
doubled  before  computing  a resultant  vector,  the  calculations  used  to 
determine  the  components  of  the  resultant  vector  are  as  follows: 


N-S  component  = Zn^  cos (20^) 


E-W  component  = Zn^  sin (20  J 


(1) 

(2) 


where:  n.  = number  of  observations  oriented  at  0. 

l l 

0i  = angle  between  particle  long  axis  and  the  N-S  direction. 

The  orientation,  0,  and  the  length,  r,  of  the  resultant  vector  may  be 
calculated: 


tan (2  0)  = 


Zn.  sin  20. 

l l 

Zn.  cos  20. 

i l 


(3) 


r = *^Zn^  sin20^) 


+ (Zn.  cos20.) 

1 l 


(4) 


The  vector  magnitude,  L,  is  defined  as  follows: 

L = X 100%  (5) 

i 

The  vector  magnitude  may  vary  from  0 to  100  percent.  A perfectly 
uniform  distribution  of  orientations  will  yield  a resultant  vector  magnitude 
of  0 percent.  A vector  magnitude  of  100  percent  indicates  that  all 
orientations  are  exactly  the  same  or,  in  the  case  of  grouped  data,  all 
orientations  lie  within  the  same  class  interval.  According  to  Curry  (1956) , 
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particles  per  thin  section.  The  values  of  L,  0 and  p were  then  calculated. 
If  the  significance  level  was  less  than  0.95  (i.e.  p > 0.05)  additional 
observations  were  made  until  the  calculated  value  of  p was  found  to  be  less 
than  or  equal  to  0.05. 

It  should  be  noted  that  a truly  uniform  distribution  of  long  axis 
orientations  does  not  lend  itself  to  this  method  of  analysis.  The 
relationship  expressed  in  Equation  8 implies  that  the  low  values  of  h which 
would  result  from  sampling  such  a uniform  distribution  would  yield  large 
values  for  p,  even  if  an  infinitely  large  number  of  observations  were 
tabulated.  As  nearly  uniform  distributions  are  seldom  encountered  in 
practice,  the  Rayleigh  test  of  significance  remains  a viable  tool  for  use 
in  sand  fabric  studies. 


Interparticle  Contact  Orientations 

The  orientation  of  any  interparticle  contact  can  be  represented  by 
1 2 

the  normals,  , to  the  tangent  plane,  as  shown  in  Fig.  15.  A method 

for  analyzing  the  orientation  of  contact  plane  normals,  developed  by  Oda 
(1972) , was  used  in  this  investigation. 

It  is  convenient  to  express  vector  orientations  in  terms  of  spherical 
coordinates.  Spherical  and  Cartesian  coordinates  are  related  as  follows: 

x = r sinB  cosd  (9a) 

y = r sin3  sina  (9b) 

z = r cosB  (9c) 


where  r is  the  vector  length  and  the  angles  a and  B describe  vector  direction 
as  shown  in  Fig.  15. 

An  element  dA  of  the  surface  area  of  a sphere  is  given  by: 

dA  a = rdB  • rsinBda  = r2  sinB  dadB  (10) 

a,  p 
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for  a sphere  of  unit  radius: 


dA^  a - sin$  dadB 

a,p 


Integration  yields  a value  for  the  entire  surface  area,  A: 

sinB  dadB 


(ID 


a 


(12) 


Consider  M contacts  within  an  assemblage  of  granular  particles.  At 

each  interparticle  contact  there  are  two  contact  surfaces,  one  belonging 

1 2 

to  each  grain.  Normal  directions  and  N^  in  Pig.  15  are  perpendicular 

2 

to  the  tangent  plane  at  the  contact  point.  The  direction  of  is 
described  by  angles  oi  and  8;  whereas,  the  direction  of  N_^  is  described  by 
the  angles  a + it  and  ir  - B. 

The  distribution  of  interparticle  contact  normals  can  be  shown  by 
plotting  their  piercing  points  on  a sphere  or  its  graphical  equivalent, 
the  equal  area  stereonet  discussed  by  Turner  and  Weiss  (1963) . The  total 
number  of  plotted  points  is  2M,  twice  the  number  of  contacts. 

The  total  number  of  contact  normals  within  the  angular  ranges  a + da 
and  B + dB  is: 


2M  E(a,B)  sinB  dadB 


(13) 


where  E(a,B)  is  the  probability  density  of  points  of  contact  within  the 
angular  intervals,  a function  of  both  a and  B.  All  plotted  points  will  be 
included  if  Equation  (13)  is  integrated  over  the  intervals  of  a = 0°  to 
360°  and  B = 0°  to  180°.  Thus  E(a,B)  must  satisfy  the  following  equality: 


// 

a 8 


E(a,B)  sin  dadB  * 1 


(14) 
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In  the  present  investigation  absolute  values  of  a are  not  known,  since 
horizontal  orientations  were  not  marked  on  the  Osterberg  samples.  Quantitative 


i 


comparisons  between  samples  have  thus  been  limited  to  consideration  of  the 
distributions  of  3^,  the  angles  between  contact  normals  and  the  vertical 
direction.  To  eliminate  the  dependance  of  E(a,$)  on  a,  E(a,B)  may  be  summed 
across  all  values  of  a.  Thus  Equation  14  becomes: 


2n 


fm, 

6 


sinBdB  = 1 


(15) 


1 2 

Since  the  vectors  JT  and  are  colinear  E(B)  equals  E (tt— 3)  . For 

example,  the  probability  density  in  the  region  B = 0°  to  3 = 10°  is  the 
same  as  in  the  region  3 = 180°  to  3 = 170° . As  a result  only  the  vector 
orientations  in  one  hemisphere  need  be  measured  and  plotted  and  the  results 
doubled  to  obtain  the  density  distribution  over  the  entire  sphere.  Thus 
Equation  15  becomes 


or 


90° 

2 JVtt  y*  E(B)  sinBdB J = 1 
0° 


sinBdB  + 


E(3)  sin3d3 


] 


(16) 
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If  is  the  total  number  of  normals  in  the  range  = 0°  to  = 10° , then: 

10° 


/ 


Mx  = 4ttm  / E(B)  sinBdB 


(18) 


Integration  gives  the  mean  value  of  E^O)  within  the  range  08  3 < 10°: 

M. 


E1(P)  = 4ttm ( cos  0°  - cos  10°) 
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(19) 
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Mean  values  of  E(3)  for  other  intervals  are  calculated  in  a similar  manner. 

For  an  isotropic  fabric  E^($)  = constant  = 1/41T.  Values  of  E^(3)  > l/4ir 
mean  a greater  than  average  proportion  of  contact  normals  at  an  inclination 
of  S..  A value  of  E^CS)  < 1/4tt  indicates  a less  than  average  proportion  of 
contact  normals  at  inclination  3^. 

Pore  Size  Measurements 

The  nature  of  pore  spaces,  as  discussed  by  Mitchell  (1976) , is  an 
important  element  of  soil  fabric.  A relatively  simple  approach  to  the 
approximation  of  pore  sizes  in  intact  samples  is  to  measure  the  equilibrium 
water  content  at  various  values  of  negative  pore  water  pressure  (soil  suction) 
The  principle  of  capillary  rise  provides  the  basis  for  inferring  pore  size 
distribution  from  measured  values  of  water  content  and  soil  suction. 

Consider  a glass  capillary  tube  of  radius,  R,  as  illustrated  in 
Fig.  16:  The  attraction  between  the  glass  and  the  water  molecules  combined 
with  the  surface  tension  of  the  water  pulls  the  water  up  into  the  tube  to 
a height,  hc>  above  the  outside  water  level.  The  height,  hc,  is  known  as 
the  "height  of  capillary  rise" . The  upper  surface  or  meniscus  of  the  water 
assumes  a concave  shape  which  joins  the  walls  of  the  tube  at  an  angle,  a, 
known  as  the  "contact  angle".  The  value  of  a depends  upon  the  material  that 
constitutes  the  capillary  wall  and  on  the  type  of  impurities  that  cover  it. 

For  quartz  tubes  with  chemically  clean  or  wetted  walls,  a is  equal  to  0°. 

If  a denotes  the  surface  tension  of  water  and  y denotes  the  unit 

'w 

weight  of  water,  equilibrium  requires  that 


whence 


h TO  Y = 2TTR0COSCI 
c w 


, 20 
h = — — cosa 
c Ry 

w 


(20) 


(21) 
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with  quartz)  and  a = 0.075  gm/cmf  Equation  (21)  may  be  rewritten  as  follows 


If  soil  voids  are  assumed  to  behave  as  a large  number  of  inter- 


an  effective  or  average  pore  radius,  R,  which  will  support  a capillary  rise 


Thus  if  a 


inferred  by  using  Equation  (22)  that  a certain  volume  of  pores  exists  within 


the  sample  with  an  effective  radius  less  than  or  equal  to  R.  If  a series  of 


the  relative  volumes  of  pore  spaces  with  different  maximum  effective  radii 


large  void  which  is  isolated  by  small  pores  will  remain  saturated  until 


pores.  The  resulting  outflow  of  water  will  be  attributed  to  a volume  of  small 


A Soiltest  A-98  Terape  Pressure  cell  was  used  in  this  study  to  apply 


flow  of  water  into  and  out  of  each  sample.  The  equilibrium  water  content  for 


each  level  of  soil  suction  was  back  calculated  from  the  initial  water  content 


VI.  APPARENT  LONG  AXIS  ORIENTATIONS 


The  distributions  of  apparent  long  axis  orientations,  as  well  as 
- 

the  direction,  0,  and  magnitude,  L,  of  the  resultant  vectors  were  deter- 
mined using  the  procedures  previously  described.  The  longitudinal  axis 
(vertical  axis)  of  the  original  cylindrical  sample  was  used  as  the  reference 
axis  for  the  vertical  sections.  An  arbitrarily  chosen  "X"  axis  was  used  as 
the  reference  axis  for  horizontal  samples.  Long  axis  orientations  are 
expressed  in  terms  of  9,  the  angle  between  the  reference  axis  and  the 
observed  particle  long  axis.  Thus,  for  a vertical  section,  a value  of 
0=0°,  180°  means  that  the  particles  are  oriented  in  a vertical  direction; 
while  a value  of  9 = 90°,  270°  indicates  a horizontal  orientation. 

The  calculated  values  of  mean  vector  orientation  and  vector  magnitude 
are  summarized  in  Table  2.  Frequency  histograms  based  on  10°  azimuth  class 
intervals  were  prepared  from  the  data  on  each  thin  section.  The  frequency 
histograms  for  each  sample  are  reproduced  in  Figs.  18  through  28. 


Long  Axis  Orientations  of  Particles  in  Horizontal  Thin  Sections 

Resultant  vector  orientations,  0,  obtained  from  the  horizontal  thin 
sections  are  seen  to  exhibit  a wide  degree  of  variation.  A complete 
interpretation  of  these  values  is  not  possible  since  horizontal  orientations 
of  the  Osterberg  samples  at  the  site  were  not  recorded  at  the  time  of  sampling. 

The  resultant  vector  magnitude,  L,  however,  is  independent  of  the  choice 
of  reference  axes,  so  a comparison  of  different  vector  magnitudes  can  provide 
information  concerning  the  relative  intensity  of  horizontal  particle 
orientation.  Resultant  vector  magnitudes  obtained  from  horizontal  sections 
are  presented  in  Fig.  29  as  functions  of  depth.  The  magnitudes  of  resultant 
vectors  obtained  from  the  River  site  are  seen  to  exhibit  a general  increase 
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with  depth,  ranging  from  values  of  approximately  10%  at  a depth  of  4.5  m.  to 
approximately  18%  at  a depth  of  13  m.  Resultant  vector  magnitudes  obtained 
from  the  Road  site  exhibit  a similar  increase  with  depth,  but  the  intensities 
of  horizontal  particle  orientation  are  seen  to  be  slightly  less  than  those 
observed  in  the  River  site  samples. 

Although  the  intensity  of  particle  orientation  in  a horizontal  plane 
has  not  yet  been  found  to  be  of  major  significance  to  liquefaction  behavior, 
the  lower  L values  exhibited  by  the  Road  site  are  in  agreement  with  results 
obtained  by  Mitchell  et  al.  (1976)  which  show  that  a more  random  distribu- 
tion of  long  axis  orientation  will  produce  higher  values  of  cyclic  strength 
in  compacted  laboratory  samples. 

Long  Axis  Orientations  of  Particles  in  Vertical  Thin  Sections 

Although  direct  comparisons  of  9 values  observed  in  the  vertical  thin 
sections  are  ambiguous  in  the  absence  of  fixed  horizontal  axes,  a qualitative 
approach  may  be  taken  to  evaluate  the  average  inclinations  of  particle  long 
axes  to  the  vertical  axis.  It  is  evident  that  preferred  orientation  of 
particles  in  the  vertical  direction  will  be  reflected  in  low  values  of  0 as 
observed  in  both  ZX  and  ZY  thin  sections.  The  average  value  of  0 observed 
in  the  two  vertical  thin  sections  may  thus  provide  a qualitative  indication 
of  this  preferred  orientation.  The  values  of  0 observed  in  the  ZX  and  ZY 
thin  sections  of  each  sample  are  presented  as  functions  of  depth  in  Fig.  30. 
In  the  estimated  zone  of  liquefaction,  the  average  values  of  0 from  the  Road 
site  samples  are  lower  than  those  of  the  River  site  samples.  Cyclic  triaxial 
strength  values  obtained  by  Silver  and  Ishihara  (1977)  are  presented  in 
Fig.  31  as  functions  of  depth.  A comparison  of  the  particle  orientations 
in  Fig.  30  with  the  strength  values  presented  in  Fig.  31  shows  that  the 
higher  cyclic  strength  of  the  Road  site  samples  for  depths  less  than  9 m 
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may  be  associated  with  a smaller  average  inclination  of  particle  long  axes 
to  the  vertical. 

In  the  absence  of  known  horizontal  reference  axes,  a quantitative 
comparison  of  vector  magnitude  (L)  values  in  different  vertical  thin 
sections  is  not  possible.  To  illustrate  this  fact,  a schematic  illustration 
of  a two  grain  sand  sample  is  presented  in  Fig.  32.  As  shown  in  the  figure, 
a thin  section  of  the  sample  prepared  in  the  ZX  plane  will  exhibit  a vector 
magnitude  which  differs  from  the  vector  magnitude  observed  in  the  thin 
section  prepared  in  the  ZY  plane.  Nonetheless,  the  average  value  observed 
in  the  ZX  and  ZY  thin  sections  may  be  used  as  a qualitative  index  of  the 
intensity  of  particle  orientation.  L values  from  vertical  thin  sections 
are  presented  as  functions  of  depth  in  Fig.  33.  The  Road  site  samples  are 
seen  to  exhibit  higher  L values  than  the  River  site  samples  implying  that 
particles  in  the  Road  site  samples  are  more  strongly  oriented.  Furthermore, 
the  variation  of  average  L values  with  depth  is  seen  to  be  similar  at  the 
two  sites. 

These  results  would  appear  to  be  contrary  to  the  results  reported  by 
Mitchell  et  al.  (1976)  which  showed  that,  for  a given  relative  density,  a 
more  random  particle  orientation  (lower  L)  produced  higher  strength  under 
cyclic  triaxial  loading.  However,  in  the  previous  study  the  conclusion 
was  drawn  from  tests  in  which  the  preferred  direction  of  long  axis 
orientations  was  in  a direction  normal  to  the  direction  of  the  applied 
cyclic  deviator  stress.  The  results  of  the  present  study  suggest  that  the 
greater  the  intensity  of  preferred  long  axis  orientation  in  a direction 
parallel  to  the  applied  cyclic  deviator  stress,  the  greater  the  resistance 
to  liquefaction.  Thus,  the  results  of  the  two  studies  may  not  be  inconsis- 
tent with  each  other,  but  reflect  the  importance  of  stress  direction  relative 
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to  orientation  direction  in  determining  resistance  to  cyclic  loading. 


The  values  of  vector  magnitude  (L) , for  long  axis  orientations  in  the 
vertical  thin  sections  (Fig.  33)  are  predominantly  in  the  range  of  10  to 
20  percent.  Recalling  that  the  average  axial  ratio  (AR)  for  the  Niigata 
sand  is  about  0.6,  reference  to  Oda  (1976)  shows  that  values  in  this  range 
are  characteristic  of  samples  prepared  by  up  and  down  vibrations.  Samples 
prepared  by  pouring,  however,  give  vector  magnitudes  considerably  larger, 
of  the  order  of  35  to  40  percent. 


VII.  INTERPARTICLE  CONTACT  ORIENTATIONS 


Interparticle  contact  orientations,  expressed  in  terms  of  the  normal 
to  the  mterparticle  contact  plane,  were  measured  in  each  thin  section  and 
plotted  on  equal  area  stereonets  according  to  the  procedure  presented  by 
Turner  and  Weiss  (1963).  The  density  function,  E($)»  was  calculated 
according  to  the  procedure  outlined  in  Section  V of  this  report. 

The  equal  area  stereonets  are  presented  in  Figs.  34  through  44.  On 
these  diagrams  blank  areas  indicate  contact  normal  directions  within  which 
less  than  one  percent  of  the  normals  intersect  one  percent  of  the  surface 
area  of  a sphere  containing  all  possible  directions.  Light  stippled  areas 
enclose  regions  containing  1 to  2 percent  of  the  contact  normals  per  one 
percent  of  the  area  of  the  circumscribing  sphere.  The  heavily  stippled 
areas  indicate  more  than  2 percent  of  the  normals  per  one  percent  of  the 
sphere  surface.  Thus  the  darker  the  area  on  the  stereonet,  the  greater  the 
concentration  of  interparticle  contact  normals  in  the  directions  represented 
by  that  area. 

It  should  be  noted  in  viewing  these  stereonets  that  because  the 
horizontal  orientation  of  the  samples  was  not  known  and  the  selection  of  X 
and  Y axis  directions  was  arbitrary,  the  stereonet  representing  the  true 
in-situ  distribution  of  interparticle  contact  normals  is  indeterminate.  In 
any  case  the  pattern  would  be  the  same  as  shown  in  the  figures  but  "rotated" 
to  the  proper  position  around  the  Z axis. 

When  the  influence  of  rotation  around  the  Z axis  is  considered,  it 
should  be  remembered  that  an  axis  through  the  center  of  a sphere  intersects 
the  surface  of  the  sphere  at  two  diammetrically  opposite  points.  If  the 
sphere  is  rotated  from  east  to  west,  an  axis  intersection  point  which 
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disappears  in  the  west,  will  immediately  have  its  diametrically  opposite 
point  of  intersection  appearing  in  the  east.  The  simultaneously  dis- 
appearing and  appearing  points  of  intersection  are  equidistant  from,  but 
on  opposite  sides,  of  the  equator. 

The  stereonets  in  Figs.  34  to  44  were  studied  by  Professor  Richard 
E.  Goodman  of  our  Geotechnical  Engineering  faculty.  Professor  Goodman  is 
an  ejqpert  on  stereonets  and  their  interpretation.  His  study  was  done  with- 
out prior  knowledge  of  where  the  samples  were  from  or  their  groupings  from 
two  different  sites.  His  conclusions  were  that  the  patterns  appear  to  fall 
into  at  least  two  groups.  One  of  them  is  represented  by  sanples  Bl-Sl 
(Fig.  34) , B1-S7  (Fig.  38)  and  B3-S5  (Fig.  43) , where  there  is  approximately 
orthogonal  synmetry  with  three  preferred  orientations.  All  the  others, 
except  possibly  those  for  sample  B1-S6  (Fig.  37)  and  sample  B3-S2  (Fig.  40) 
indicate  non-orthogonal  synmetry  with  more  than  three  preferred  orientations. 

An  important  conclusion  that  results  from  the  interpretations  by 
Professor  Goodman  is  that  the  stereonets  do  not  indicate  significant 
difference  between  the  interparticle  contact  normal  orientations  for  River 
site  (B1  samples)  and  the  Road  site  (B3  samples) . 

This  lack  of  significant  differences  between  the  distributions  of 
contact  normals  for  samples  from  the  two  sites  is  borne  out  also  by  the  dis- 
tribution functions  E(0)  shown  in  Figs.  45  and  46.  The  distributions  do 
suggest  that  the  concentration  of  interparticle  contact  normals  in  a near 
vertical  direction  (contact  planes  in  a horizontal  direction)  increases  with 
depth  for  both  sites.  Oda  (1976)  has  shown  that  an  increase  in  the  concen- 
tration of  normals  in  the  maximum  principal  stress  direction  is  related  to 
an  increase  in  strength  of  a granular  material.  This  is  because  contact  planes 
perpendicular  to  the  major  stress  direction  are  the  most  effective  to  carry 


that  stress. 
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VIII.  SOIL  SUCTION-WATER  CONTENT  RELATIONSHIPS 


As  described  in  Section  III,  portions  of  the  frozen  samples  were  cut, 
trimmed  and  allowed  to  thaw  in  the  Tempe  Pressure  cell.  After  a thawing 
period  of  4 hours  successive  increments  of  negative  pore  water  pressure  (soil 
suction)  ranging  from  1 cm.  to  900  cm.  of  water  head  were  induced  and  the 
amount  of  water  which  flowed  from  the  sample  was  measured.  The  first  phase 
of  the  test,  in  which  increasing  values  of  soil  suction  were  applied,  is 
termed  the  "drying  cycle".  The  second  phase  of  the  test,  in  which  the  soil 
suction  was  decreased  in  successive  increments,  is  termed  the  "wetting  cycle." 
The  equilibrium  sample  water  content  corresponding  to  each  value  of  applied 
soil  suction  could  be  calculated  from  the  initial  and  final  sample  weights 
and  the  amount  of  water  which  flowed  in  or  out  of  the  sample. 

Sample  water  contents  are  plotted  against  the  corresponding  soil 
suction  values  in  Figs.  47  through  57.  Initial  sample  water  contents,  are 
seen  to  fall  within  the  range  of  8 to  10  percent  for  all  samples.  The  results 
show  that  most  of  the  water  outflow  measured  during  the  drying  cycles  occurred 
in  the  soil  suction  range  of  10  to  60  cm.  H20. 

The  results  for  all  the  River  site  samples  have  been  averaged  and  are 
shown  in  Fig.  58.  An  average  water  content  vs.  suction  curve  for  the  Road 
site  samples  is  shown  in  Fig.  59.  The  behavior  of  the  samples  from  the  two 
sites  is  compared  in  Fig.  60  in  terms  of  degree  of  saturation  vs.  soil 
suction  and  vs.  equivalent  pore  radius  as  calculated  from  equation  (22). 

It  is  clear  from  Fig.  60  that  the  River  site  samples  had  a higher 
degree  of  saturation  at  a given  value  of  soil  water  suction  than  did  the 
Road  site  samples.  As  the  tests  followed  a drying  cycle  (increasing  suction) 
followed  by  a wetting  cycle,  this  implies  that  a greater  proportion  of  the 


void  space  was  made  up  of  small  pores  in  the  River  site  samples  than  in  the 
Road  site  samples.  A part  of  this  difference  may  have  resulted  from  the 
fact  that  the  Road  site  samples  were  less  dense  (average  void  ratio  of  0.91) 
than  the  River  site  samples  (average  void  ratio  of  0.84)  used  for  the  suction 
tests. 

Nonetheless,  the  degree  of  saturation  represents  the  proportion  of 
the  void  space  that  is  filled  with  water,  and  the  ordinate  values  in  Fig.  60 
can  be  interpreted  as  percentage  of  pores  finer  than  any  given  size.  It 
is  clear  that  the  River  site  samples  contain  more  small  pores  than  do  the 
Road  site  sanples.  This  may  be  a factor  in  the  apparent  lower  resistance  to 
liquefaction  at  the  River  site  than  at  the  Road  site  in  spite  of  the  higher 
void  ratio  of  the  Road  site  samples.  A high  proportion  of  small  pores  means 
a low  hydraulic  conductivity.  A low  hydraulic  conductivity  means  limited 
ability  for  pore  pressures  generated  during  cyclic  loading  to  dissipate 
rapidly,  and,  hence,  a potentially  more  rapid  onset  of  liquefaction. 
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IX.  X- RADIOGRAPHS 


As  described  in  Section  III,  remaining  portions  of  the  resin  impreg- 
nated samples,  after  thin  section  preparation , were  prepared  for  X- 
radiographic  inspection.  Actual-size,  positive  X-radiographs  of  horizontal 
and  vertical  sections  cure  shown  for  samples  B1-S3,  B1-S4,  B1-S5,  B1-S6  and 
B1-S7  (all  from  the  river  site)  in  Fig.  61,  while  horizontal  and  vertical 
sections  for  samples  B3-S1,  B3-S2,  B3-S4,  B3-S5  and  B3-S7  (all  from  the 
road  site)  cure  shown  in  Fig.  62. 

Comparison  of  Figs.  61  and  62  shows  that  the  samples  taken  at  9.5  m 
depth  (B1-S5  and  B3-S5)  appear  reasonably  similar  as  far  as  grain  size  is 
concerned.  The  same  holds  true  for  the  samples  taken  at  12.5  m (B1-S7  and 
B3-S7)  depth,  while  the  river  site  sample  taken  at  4.5  m (B1-S3)  contains 
larger  particles  than  the  road  site  sample  from  a similar  depth  (B3-S1) . 

The  larger  grains  in  sample  B1-S3,  when  compared  to  B3-S1,  are  not  reflected 
in  the  grain  size  distribution  curves  for  these  samples  (Figs.  6 and  7) . 

A slightly  darker  zone  can  be  observed  on  the  circular  edges  of  all 
horizontal  sections.  The  darker  zone  is  probably  caused  by  densification 
during  sampling  of  a narrow  annulus  of  soil  immediately  adjacent  to  the  sampler. 
The  disturbed  zone  can  be  seen  to  be  approximately  1/8  to  1/4  inch  (3  to 
6 mm)  thick. 

The  variation  in  X-radiograph  intensity  shown  by  samples  B1-S3, 

B1-S5,  B3-S1  and  B3-S5  may  be  indicative  of  variations  in  density  of  these 
samples  on  the  scale  of  the  sections.  It  would  not  be  possible  to  quantita- 
tively assess  these  variations  by  means  of  relative  density  measurements, 
as  the  latter  require  relatively  large  samples. 
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Similarly  a zone  of  lower  density,  about  1/8  inch  (3  mm)  thick,  is 


i 

also  shown  by  the  vertical  section  of  sample  B1-S4. 

Sanple  B1-S6,  from  a depth  of  11.5  m,  was  located  near  the  assumed 
maximum  depth  of  liquefaction  (9.0  m) . The  striations  in  the  vertical 
section  may  have  been  caused  by  shear  displacements. 

The  cracks  shown  by  light  zones  in  the  horizontal  sections  of  samples 
B3-S1  and  B3-S4  are  probably  not  representative  of  in-situ  conditions,  and 
may  have  been  caused  by  either  the  sampling  or  resin  impregnation  processes, 
while  the  darker  lines  in  the  horizontal  section  of  sample  B3-S4  was  caused 
by  sawing  as  evidenced  by  inspection  of  the  section  surface. 

While  the  X-radiographs  provide  a reasonably  good  visual  representa- 
tion of  the  respective  samples,  the  number  and  size  of  samples  is  insuf- 
ficient to  allow  any  conclusion  regarding  the  observed  different  lique- 
faction behavior  at  the  two  sites. 
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X.  SUMMARY  AND  CONCLUSIONS 


The  fabric  of  11  undisturbed  sand  samples  from  Niigata,  Japan, 
obtained  using  an  Osterberg  sampler,  has  been  studied  using  petrographic 
thin  sections  and  water  content-soil  water  suction  measurements.  Six 
samples,  from  a depth  range  of  4.5  to  13.5  m,  were  obtained  from  a single 
boring  (Bl)  at  a location  termed  the  River  site,  where  the  sand  is  known 
to  have  liquefied  during  the  1964  Niigata  earthquake.  Five  samples  from 
a depth  range  of  4.5  to  13.5  m were  obtained  from  a single  boring  (B3)  at 
a location  termed  the  Road  site,  where  no  liquefaction  was  observed  during 
the  1964  earthquake.  All  samples  were  frozen  after  they  were  taken  from 
the  ground  prior  to  shipment  to  the  University  of  California,  Berkeley, 
Geotechnical  Engineering  Laboratory  for  study. 

The  major  findings  from  this  investigation  are  as  follows. 

1.  The  techniques  developed  for  preparation  of  thin  sections  for 
fabric  study  and  for  samples  for  soil  water  content  vs.  suction 
measurements  worked  well. 

2.  The  gradational  and  particle  shape  characteristics  of  the  medium 
to  fine  sands  at  the  Road  and  River  sites  are  similar,  although 
the  average  particle  size,  D5Q,  and  uniformity  coefficient, 
D60/Dl0»  at  the  R*ver  site  are  slightly  larger  than  at  the  Road 
site.  Particle  shape,  described  in  terms  of  the  average  particle 
axial  ratio  (AR)  is  quite  similar  for  the  two  sites;  0.62  for 
the  River  site  and  0.59  for  the  Road  site,  with  the  River  site 
sanples  exhibiting  a slightly  wider  range  of  values. 
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3.  The  void  ratios  of  the  samples  are  distributed  erratically 
with  depth  at  both  sites;  however,  the  void  ratios  of  the  Road 
site  samples  are  generally  higher  than  those  from  the  River 
site. 

4.  Preferred  orientation  of  particle  long  axes  in  any  known 
horizontal  direction  could  not  be  determined  owing  to  the  absence 
of  sample  orientation  information.  The  intensity  of  horizontal 
orientation  relative  to  arbitrary  axes  was  determined,  however, 
in  terms  of  the  resultant  vector  magnitude,  L.  For  samples  from 
the  upper  9 m,  the  presumed  zone  of  liquefaction  at  the  River  site, 
the  Road  site  samples  exhibited  slightly  lower  L values  (Fig.  29) , 
indicative  of  more  random  particle  arrangement.  This  is  consistent 
with  the  results  of  a previous  study  (Mitchell  et  al.,  1976)  which 
indicated  that  samples  with  a more  random  orientation  of  particles 
in  the  horizontal  plane  were  more  resistant  to  liquefaction  in 
cyclic  triaxial  tests. 

5.  The  mean  orientation  direction  9 of  particle  long  axes  in  vertical 
thin  sections  was  closer  to  the  vertical  for  samples  from  the  upper 

9 m from  the  Road  site  than  for  samples  from  the  River  site  (Fig.  30) . 

6.  The  magnitudes  of  L observed  in  different  vertical  thin  sections 
provide  qualitative  indices  of  the  degree  of  particle  orientations 
in  vertical  planes.  The  Road  site  samples  exhibit  a stronger 
degree  of  particle  long  axis  orientation  in  vertical  planes  than 
do  the  River  site  sanp>les.  This  finding  is  shown  to  be  consistent 
with  a higher  resistance  to  liquefaction  in  cyclic  triaxial  tests 
for  the  Road  site  samples. 
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7.  The  values  of  vector  magnitude  L for  long  axis  orientations  in 
the  vertical  thin  sections  are  characteristic  of  reconstituted 
samples  of  similar  sands  prepared  by  up  and  down  vibrations  (Oda, 
1976)  but  not  of  samples  prepared  by  pouring. 

8.  Analysis  of  equal  area  stereonets  does  not  indicate  significant 
differences  between  the  orientations  of  normals  to  interparticle 
contact  planes  for  Road  and  River  site  samples.  This  lack  of 
significant  differences  in  the  distributions  of  interparticle 
contact  normals  for  samples  from  the  two  sites  is  borne  out  also 
by  the  E(8)  distribution  functions.  At  both  sites,  however,  the 
concentration  of  normals  to  interparticle  contacts  in  a near 
vertical  direction  increases  with  depth. 

9.  The  River  site  samples  had  a higher  degree  of  saturation  at  a 
given  value  of  soil  water  suction  than  did  the  Road  site  samples. 
These  results  suggest  that  the  River  site  samples  had  a greater 
volume  of  small  size  pores  than  did  the  Road  site  samples.  This 
could  be  a factor  in  the  lower  resistance  to  liquefaction 
exhibited  by  the  River  site  samples. 

10.  X-radiography  of  slices  of  resin-impregnated  samples  would  appear 
to  offer  good  potential  for  studying  fabric  characteristics  on  a 
scale  of  a few  mm. 

In  interpreting  these  findings  it  should  be  noted  firstly  that  the 
samples  were  taken  from  the  ground  before  freezing.  Any  disturbance  that 
may  have  occurred  during  this  sampling  will  be  reflected  in  the  measured 
fabric  characteristics.  As  the  magnitude  of  this  disturbance  is  unknown, 
there  is  uncertainty  as  to  whether  the  observed  fabric  accurately  reflects 


the  true  in-situ  fabric. 
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Secondly,  the  horizontal  orientation  of  the  samples  was  not  recorded. 
Consequently  the  possible  existence  of  any  macroscopic  horizontal  preferred 
orientations  of  particles  or  particle  contacts  at  the  two  sites  cannot  be 
established. 

Thirdly,  the  River  site  samples  came  from  a zone  which  had  liquefied 
during  the  1964  earthquake.  Thus  present  difference  between  the  sand  fabric 
characteristics  at  the  two  sites  cannot  be  used  to  explain  the  original 
differences  in  cyclic  shear  strength  and  resistance  to  liquefaction. 

Finally,  the  fact  that  liquefaction  occurred  at  the  River  site  but 
not  at  the  Road  site  during  the  1964  Niigata  earthquake  could  be  due  to 
causes  other  than  differences  in  sand  fabric  or  density.  Specifically,  it 
is  noted  that  the  water  table  at  the  Road  site  is  at  a depth  of  2.5  m, 
whereas  it  is  at  0.6  m at  the  River  site,  a fact  that  would  give  a greater 
resistance  to  liquefaction,  for  a given  intensity  of  shaking,  at  the  Road 
site.  Furthermore,  whether  the  intensities  of  shaking  at  the  two  sites 
were  the  same  is  unknown. 

Consequently,  the  results  of  this  study,  taken  by  themselves,  cannot 
be  used  to  establish  conclusively  the  specific  role  of  sand  fabric  in-situ 
in  determining  resistance  to  liquefaction  or  the  changes  in  fabric  that 
acconq>any  liquefaction.  They  do  suggest,  however,  that  fabric  in-situ  is 
related  to  susceptibility  to  liquefaction  as  borne  out  by  the  differences 
in  fabric  observed  for  the  two  sites  and  the  corresponding  differences  in 
cyclic  triaxial  strength  reported  by  Silver  and  Ishihara  (1977) . Further- 
more, the  fabric  relationships  observed  herein  and  corresponding  cyclic 
strength  behavior  reported  by  Silver  and  Ishihara  appear  consistent  with 
previous  fabric  property  correlations  reported  by  Oda  (1976)  and  Mitchell 
et  al.  (1976) . 
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Further  studies  are  indicated  to  evaluate  the  effect  of  sampling  on 
fabric,  to  compare  the  in-situ  fabric  of  the  Niigata  sands  with  that  of 
reconstituted  samples  of  the  same  sand,  and  to  integrate  the  results  of 
the  present  investigation  more  completely  with  the  results  of  fabric  and 
liquefaction  studies  on  other  sands. 
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TABLE  1 


FLOWCHART  OF  RESEARCH  ON  UNDISTURBED  SAMPLING 
AND  TESTING  OF  SANDS  FROM  NIIGATA,  JAPAN 
(Adapted  from  Silver  and  Ishihara,  1977  ) 


UNDISTURBED  SOIL  SAMPLING 

1.  Large  diameter  sampler 
(Japanese) 

2.  Osterberg  sampler  (U.S.) 


FIELD  HANDLING  TECHNIQUES 

1.  Sample  drainage 

2.  Preparation  of  small  specimens 

3.  Field  freezing 


FABRIC  SAMPLES 

1.  Specimen  thawing 

2.  Epoxy  replacement 

3.  Density  and  fabric  study 

4.  Soil-water  suction 
measurement 


1 

1 

LABORATORY  TESTS  ON 

1 

1 

1 

RECONSTITUTED  SPECIMENS 

1 

1 

.1.  Study  of  effect  of 

1 

1 

sample  preparation 

1 

1 

I 

technique 

1 

COMPARISON  OF  N VALUES, 
IN-SITU  DENSITY  AND 
RELATIVE  DENSITY  WITH 
CYCLIC  SOIL  STRENGTH 


■ studies  described  in  this  report 
indicates  work  not  yet  completed 

indicates  work  described  by  Silver  and  Ishihara  (1977) 


FIG.  2 SOIL  PROFILES  FOR  THE  ROAD  AND  RIVER  SITES 
(After  Silver  and  Ishihara,  1977) 
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FIG.  4 THIN  SECTION  LOCATIONS  FROM  WITHIN  RESIN 
IMPREGNATED  SAMPLES 
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FIG.  7 GRAIN  SIZE  DISTRIBUTIONS  OF  ROAD  SITE  SAMPLES 
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FIG.  8 MEAN  PARTICLE  SIZE,  D50>  AND  UNIFORMITY  COEFFICIENT,  D^-./D.  AS  FUNCTIONS  OF  DEPTH 
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FIG.  9 THE  MEASUREMENT  OF  PARTICLE  AXIAL  RATIO 
(After  Oda,  1976) 
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FIG.  10b  PARTICLE  AXIAL  RATIO  DISTRIBUTIONS  FOR  ROAD 
SITE  SAMPLES 


FIG.  11  AVERAGE  AXIAL  RATIO  DISTRIBUTIONS  FOR 
RIVER  AND  ROAD  SITE  SAMPLES 
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o— □ Range  of  values  for  Osterberg  and 
large  diameter  samples 
(Silver  ond  Ishihara,  1977) 

o From  samples  used  for  thin  section 
preparation 

A From  samples  used  for  soil  suction  test 


FIG.  12  RIVER  AND  ROAD  SITE  SAMPLE  VOID  RATIOS  AS  FUNCTIONS 
OF  DEPTH 


IN- SITCJ  RELATIVE  DENSITY  VALUES  FOR  THE  ROAD  AND  RIVER 
SITES  AS  A FUNCTION  OF  DEPTH 
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FIG.  20  HISTOGRAMS  OF  APPARENT  PARTICLE  LONG  AXIS  ORIENTATIONS  FOR 
RIVER  SITE  SAMPLE  B1-S5 
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FIG.  23  HISTOGRAMS  OF  APPARENT  PARTICLE  LONG  AXIS  ORIENTATIONS 
FOR  ROAD  SITE  SAMPLE  B3-S1 
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FIG.  24  HISTOGRAMS  OF  APPARENT  PARTICLE  LONG  AXIS  ORIENTATIONS 
FOR  ROAD  SITE  SAMPLE  B3-S2 


. 25  HISTOGRAMS  OP  APPARENT  PARTICLE  LONG  AXIS  ORIENTATIONS 
FOR  ROAD  SITE  SAMPLE  B3-S3 
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FIG.  26  HISTOGRAMS  OF  APPARENT  PARTICLE  LONG  AXIS  ORIENTATIONS 
FOR  ROAD  SITE  SAMPLE  B3-S4 
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FIG.  27  HISTOGRAMS  OF  APPARENT  PARTICLE  LONG  AXIS  ORIENTATIONS 
FOR  ROAD  SITE  SAMPLE  B3-S5 


% Frequency 


f : RESULTANT  VECTOR  MAGNITUDES  OBTAINED 

'ROM  HORIZONTAL  THIN  SECTIONS  AS  A 
M NOTION  OF  DEPTH 


Stress  Ratio  to  Couse  Failure  in  20  Cycles 
0 0.1  0.2  0.3  0.4 


FIG.  31  COMPARISON  OF  IN-SITU  CYCLIC  TRIAXIAL 
STRENGTH  (FAILURE  IN  20  CYCLES)  AS  A 
FUNCTION  OF  DEPTH  AT  THE  ROAD  AND 
RIVER  SITES 

(From  Silver  and  Ishihara,  1977) 
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FIG.  32  ILLUSTRATION  OF  THE  DEPENDENCE  OF  OBSERVED  VECTOR 
MAGNITUDES,  L,  ON  THIN  SECTION  ORIENTATION 
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FIG.  47  EQUILIBRIUM  WATER  CONTENT  VS.  SOIL  SUCTION  FOR  RIVER  SITE  SAMPLE 
B1-S3 
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FIG.  48  EQUILIBRIUM  WATER  CONTENT  VS.  SOIL  SUCTION  FOR  RIVER  SITE  SAMPLE 
B1-S4 
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FIG.  51  EQUILIBRIUM  WATER  CONTENT  VS.  SOIL  SUCTION  FOR  RIVER  SITE  SAMPLE 
B1-S7 


Road  Site 


FIG.  52  EQUILIBRIUM  WATER  CONTENT  VS.  SOIL  SUCTION  FOR  ROAD  SITE  SAMPLE 
B3-S1 


Rood  Site 


EQUILIBRIUM  WATER  CONTENT  VS.  SOIL  SUCTION  FOR  ROAD  SITE  SAMPLE 
B3-S2 


Rood  Site 


0.25 


FIG.  57  EQUILIBRIUM  WATER  CONTENT  VS.  SOIL  SUCTION  FOR  ROAD  SITE  SAMPLE 
B3-S7 
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FIG.  59  AVERAGE  EQUILIBRIUM  WATER  CONTENT  VS.  SOIL  SUCTION  FOR 
ROAD  SITE  SAMPLES 
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